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Abstract 

Continuing work initiated in earlier publications [Yamada, Asada, Phys. Rev. D 82, 104019 
(2010), 83, 024040 (2011)], we investigate the post-Newtonian effects on Lagrange's equilateral tri- 
angular solution for the three-body problem. For three finite masses, it is found that the equilateral 
triangular configuration satisfies the post-Newtonian equation of motion in general relativity, if and 
only if all three masses are equal. When a test mass is included, the equilateral configuration is 
possible for two cases: (1) one mass is finite and the other two are zero, or (2) two of the masses 
are finite and equal, and the third one is zero, namely a symmetric binary with a test mass. The 
angular velocity of the post-Newtonian equilateral triangular configuration is always smaller than 
the Newtonian one, provided that the masses and the side length are the same. 

PACS numbers: 04.25.Nx, 95.10.Ce, 95.30.Sf, 45.50.Pk 
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I. INTRODUCTION 



The three-body problem in the Newton gravity represents classical problems in astronomy 
and physics (e.g, llH3|). In 1765, Euler found a collinear solution for the restricted three- 
body problem that assumes one of three bodies is a test mass. Soon after, his solution 
was extended for a general three-body problem by Lagrange, who also found an equilateral 
triangle solution in 1772. Now, the solutions for the restricted three-body problem are called 
Lagrange points L±, L2, L3, L4 and L5, which are well known and described in textbooks of 
classical mechanics Q. 

Lagrange points have recently attracted renewed interests for relativistic astrophysics, 
where they have discussed the gravitational radiation reaction on L 4 and L 5 analytically J4] 
and by numerical methods (5), 0] . 

As a pioneering work, Nordtvedt pointed out that the location of the triangular points is 
very sensitive to the ratio of the gravitational mass to the inertial one 7|. Along this course, 
it is interesting as a gravity experiment to discuss the three-body coupling terms at the post- 
Newtonian order, because some of the terms are proportional to a product of three masses 
as Mi x M 2 x M 3 . Such a triple product can appear only for relativistic three (or more) body 
systems but cannot for a relativistic compact binary nor a Newtonian three-body system. 

The relativistic perihelion advance of the Mercury is detected only after much larger 
shifts due to Newtonian perturbations by other planets such as the Venus and Jupiter are 
taken into account in the astrometric data analysis. In this sense, effects by the three- 
body coupling are worthy to investigate. Nevertheless, most of post-Newtonian works have 
focused on either compact binaries because of our interest in gravitational waves astronomy 
or N-body equation of motion (and coordinate systems) in the weak field such as the solar 
system (e.g. [8|). Actually, future space astrometric missions such as Gaia require a 

general relativistic modeling of the solar system within the accuracy of a micro arc-second 
Furthermore, a binary plus a third body have been discussed also for perturbations of 
gravitational waves induced by the third body 12H15]. 

The theory of general relativity is currently the most successful gravitational theory 
describing the nature of space and time. Hence it is important to take account of general 
relativistic effects on three-body configurations. The figure-eight configuration that was 
found decades ago [if], Q has been numerically studied at the first post-Newtonian [.lgj and 
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also the second post-Newtonian orders |19j . According to their numerical investigations, the 
solution remains true with a slight change in the figure-eight shape because of relativistic 
effects. 

On the other hand, the post-Newtonian collinear configuration has been recently obtained 
as a relativistic extension of Euler's collinear one, where three bodies move around the 



common center of mass with the same orbital period and always line up [20|. It may offer 
a useful toy model for relativistic three-body interactions, because it is tractable by hand 
without numerical simulations. The uniqueness of the collinear configuration has been also 



proven 



2l|. 



Lagrange's equilateral triangular solution has also a practical importance, since it is stable 
for some cases. Lagrange's points L 4 and L 5 for the Sun- Jupiter system are stable and indeed 
the Trojan asteroids are located there. Clearly it is of greater importance to investigate 
Lagrange's equilateral triangular solution in the framework of general relativity. Do the post- 
Newtonian effects admit such a triangular solution? No one doubts whether the particular 
configuration is still possible at the post-Newtonian order. We shall study this issue in this 
paper. The main purpose of this paper is to show that the equilateral triangular configuration 
can satisfy the post-Newtonian equation of motion, if and only if all three finite masses are 
equal. Throughout this paper, we take the units of G — c = 1. 

II. NEWTONIAN LAGRANGE'S EQUILATERAL TRIANGULAR SOLUTION 

First, we consider the Newton gravity among three masses denoted as Mj (I = 1,2,3). 
The location of each mass is written as Xj. We choose the origin of the coordinates, so that 

M x xx + M 2 x 2 + M 3 x 3 = 0. (1) 

We start by seeing whether the Newtonian equation of motion for each body can be satisfied 
if the configuration is an equilateral triangle. Let us put R\ 2 = R 23 = R31 = a, where we 
define the relative position between masses as 

R I} = xj - x J: (2) 

and Rjj = \Ru\ for I, J — 1, 2, 3. Then, the equation of motion for each mass becomes 

d2Xl M% (3) 



alt 2 a 3 
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where M denotes the total mass ^ 7 Mj. Therefore, it is possible that each body moves 
around the common center of mass with the same orbital period. Eq. fl3]) gives 

M 



N 



(4) 



where denotes the Newtonian angular velocity. 

Figure [TJ shows an equilateral triangular configuration. Let tj denote the relative position 
vector of each mass with respect to the common center of mass (but not the geometrical 
center of the triangle) in the corotating frame with the angular velocity u^. The angles 
between (£ 2 ,£i), (l3,£ 2 ) and (£1,^3) are defined cyclically as 9i, 9 2 and 6*3 respectively. 
They are constant with time. There is an identity as 9\ + 82 + 6*3 = 2tt, which can be used 
to delete one of the angles. The orbital radius £1 = \£j\ of each body with respect to the 
common center of mass is obtained as {2] 



— \/M 2 2 + M 2 M 3 



a 

Tl 

a 

M 
a 



Ml, 



ij\lMl + M 1 M 3 + Mi, 



M 



Ml + M 1 M 2 + Ml 



(5) 
(6) 
(7) 



III. POST-NEWTONIAN EQUILATERAL TRIANGULAR SOLUTION 



Next, we consider the post-Newtonian effects on the triangu 



the Einstein- Infeld- Hoffman (EIH) equation of motion as 22H24] 



ar configuration by employing 
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where vj denotes the velocity of each mass in an inertial frame and we define 
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(8) 
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Let us see whether the three masses at the apices of an equilateral triangle can satisfy 
the EIH equation of motion. For such an equilateral triangle case, the second-order-mass 
terms are easy to handle, because every Rjj is the same as a. What we have to take care of 
is the velocity-dependent terms. 

We consider three masses in circular motion with the angular velocity co, so that each £j 
can be a constant. The position and velocity of each body are expressed as 



a?i = £i\ , (10) 




vi = £iu\ , (11) 



. cos (cot + 9i) , 

x 2 = £2 I 1 I , (12) 

sin (cot + 61) 

, — sin (cot + 61) , 
v 2 = £ 2 co\ 1 I , (13) 

cos (cot + 61) 

, cos (cot — 63) . 
xs = £3 I I , (14) 

sin (cot — 9 3 ) 

, — sin (cot — 63 1 , 
v 3 = £ 3 co I 1 I , (15) 

cos (cot — 9 3 ) 

where we used 9\ + 9 2 = 2n — 9 3 , \xi\ = £1 and \vj\ = £jco. 

For the later convenience, we compute the inner products between the velocity and rela- 
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tive position vectors as 



Rl2 


■ Vx = 


—£i£ 2 u sin 9i, 






(16) 


R12 ■ 


' v 2 = 


—t\t 2 uj sin 9i, 






(17) 


R12 ■ 


■v 3 = 


t 3 t\UJ sin #3 — 


£ 2 £ 3 co sin ( 




(18) 


R23 ■ 


■ v x = 


t\t 2 uj sin 61 — 


£3^10; sin ( 


0i + 2 ), 


(19) 


R23 ■ 


■ v 2 = 


-£ 2 £ 3 oo sin 2 , 






(20) 


R23 ■ 


■ v 3 = 


-£ 2 £ 3 u sin 9 2l 






(21) 


R31 ■ 


■ Vx = 


—£ 3 £iu sin #3, 






(22) 


R31 ' 


■ v 2 = 


£ 2 £ 3 ui sin 9 2 — 


£\£ 2 uj sin ( 


^2 + ^3), 


(23) 


R31 ■ 


W 3 = 


—£ 3 £\uj sin 3 , 






(24) 


v 1 ■ 


■v 2 = 


£i£ 2 u 2 cos 0i, 






(25) 


v 2 ■ 


■ v 3 = 


£ 2 £ 3 oo 2 cos 2 , 






(26) 


V 3 ■ 


■ Vx = 


£ 3 £\uj 2 cos 3 . 






(27) 



Note that the inner product is defined in the 3-dimensional Euclid space as ex. ■ (3 = Sija 1 ^ 
for i,j = 1,2,3. This is because the terms expressed by Eqs. ffT6"j) - fl2~T|) appear at the post- 
Newtonian order but not at the Newtonian one. Therefore, it is sufficient to use the flat-space 
metric for computing these post-Newtonian terms. Corrections to the inner product in a 
curved spacetime make 2PN (or higher order) contributions and they are thus ignored in 
this paper. 

Let us explain how to describe an equilateral triangular configuration in the Newtonian 
frame and the relativistic one. The same side length is given for the Newtonian triangle 
and the post-Newtonian one, so that each mass position can be denoted as the same vector 
Xi for both cases, whereas the position of mass center (not a geometric center) may be 
different. This treatment makes computations easier than another way that assumes the 
mass center position to be fixed as the same location for both frames by introducing three 
post-Newtonian position vectors xf N (for the three masses) possibly different from the three 
Newtonian position vectors xf . Three position vectors are needed to calculate in the latter 
approach, whereas in this paper one position vector denoting the mass center is sufficient to 
specify the system. 

In order to compute the orbital radius of each mass, the location of the mass center at 
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the post-Newtonian order must be determined. It is expressed as [22l. |23| 
G PN = E~ l J2 m a*a 

A 

where E is defined as 



1 / 2 r M B 



(28) 



i + ; vi 



M, 



(29) 



By using Eq. f[28|) . the post-Newtonian orbital radius £pjvi 



PN\ 



$ij( x \ — G l PN )(x{ — G J PN ) is obtained as 



-PNl 
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+ M 2 M| - 2M\M 2 M 3 + M x MlM z + M X M 2 M. 



3, ,2 



1 



a Co> 



M 



(30) 



By noting Eq. (JIJ), we find that the second term in the R.H.S. of Eq. ( |30l) vanishes and 
hence £pm — 1\- By the cyclic permutations, we find also tpNi = ^2 and £pN3 = ^3- 

As a consequence, the common center of mass for the equilateral solution remains un- 
changed. Without this unexpected thing, our calculations would become much more lengthy. 
The above expressions for the inner products are substituted into the R.H.S. of Eq. flE]). 
After straightforward calculations, the equation of motion for M\ can be written as 

M 



CU 2 Xi 



-x 1 + gp N iXi 



ar 
16a 3 



M 2 M 3 {M 2 -M 3 ) 
' n±1 Mf + M 2 M 3 + Ml 
a 3 



10 + ^ 2 -(-4M 1 + 5M 2 + 5M 3 ) 



(31) 



where we used Eq. (Tj0) for velocity-dependent terms, n±i = V\ji\bj is defined as the unit 
normal vector to X\, and gpNi denotes the post-Newtonian terms defined as 

1 



gpNi 



16a 4 (Mf + M 2 M 3 + Mf 

X 



80(M 2 2 + M 2 M 3 + M 3 2 )M 2 - 2(8M 2 3 + M 2 2 M 3 + M 2 M 3 2 + 8M 3 )M 
32(M 2 + M 2 M 3 + M 3 2 )M 2 + Y1M 2 M 3 {M 2 + M 3 )M 

-(16M 2 4 + 41M 3 M 3 + 84M 2 M 3 2 + 41M 2 M 3 + 16M 3 4 )}^^ 



(32) 



Here, terms with u% come from the velocity-dependent terms and may be reexpressed by 
using Eq. (TJ]). 

We should note that the third term in the R.H.S. of Eq. (131 p is parallel to the velocity of 
Mi and thus perpendicular to Xi for a circular motion case. Therefore, the mass Mi can be 
in circular motion, if and only if the coefficient of the third term vanishes, that is M 2 = M 3 . 
Likewise, the masses Mi and M 3 can be in circular motion, if and only if M 3 = Mi and 
Mi = M 2 , respectively. Hence, all the three masses can have a circular motion, if and only 
if Mi = M 2 = M 3 . 

The previous paragraph postulates that three masses are finite. Here, a test mass is 
included. First, we consider a case that two of the masses are finite and one is zero, for 
instance Mi — > 0. The third term in the R.H.S. of Eq. (1311) vanishes, if and only if M 2 = M 3 . 
Clearly, the corresponding terms of the equation of motion for M 2 and M3 vanish in the 
limit as Mi — > 0. Hence, the equilateral configuration is possible for an equal-mass binary 
with a test mass. Next, we consider a case that one mass is finite and the other two are zero, 
for instance M 2 — >■ and M 3 — > 0. When we put 0(e) = 0(M 2 ) = 0(M 3 ), the third term 
in the R.H.S. of Eq. ( 13~TT) becomes 0(e) and thus vanishes as e — > 0. The corresponding 
terms for M 2 and M 3 also become 0(e) and thus vanish as e — > 0. Therefore, the equilateral 
configuration is possible also for one finite mass and two test masses. This case is reasonable, 
since it is nothing but two test masses orbiting around the Schwarzschild black hole (in the 
weak-field approximation). 

The remaining thing to do is to see whether orbital periods of the three masses are all 
the same in order to preserve the triangular shape if Mi = M 2 = M 3 . It is easy to see 
this, because one can obtain the post-Newtonian forces gpm and gpN3 from Qpni by cyclic 
manipulations asl— > 2 — > 3 — > 1, and finally by taking the equality of Mi = M 2 = M 3 , 
one can find gpm = 9pn2 = 9pnz- Therefore, it is concluded that the equilateral triangular 
configuration remains true for the post-Newtonian equation of motion in general relativity, 
if and only if all three masses are equal. 

Eq. ( 13T|) gives uniquely the post-Newtonian angular velocity as u 2 = Ma~ 3 — gpN, where 
9pn = 9pni = 9pn2 = gpm for Mi = M 2 = M 3 = M/3. Here, g PN simply becomes 



One can show gpN > and hence 00 < oon- This means that the angular velocity of the 




(33) 



8 



M 2 




FIG. 1: Equilateral triangular configuration. Each mass is located at one of the apices. We define 
8 j (I = 1, 2, 3) with respect to the common center of mass. 

post-Newtonian equilateral triangular configuration is always smaller than the Newtonian 
one, provided that the masses and the side length are the same. This behavior occurs also 
in the post-Newtonian collinear configuration [2l|. 

IV. SUMMARY 

We investigated the post-Newtonian effects on Lagrange's equilateral triangular solution 
for the three-body problem. For three finite masses, we found that the equilateral triangular 
configuration satisfies the post-Newtonian equation of motion in general relativity, if and 
only if all three masses are equal. When a test mass is included, the equilateral configuration 
is possible for two cases: (1) one mass is finite and the other two are zero, or (2) two of the 
masses are finite and equal, and the third one is zero, namely a symmetric binary with a 
test mass. 

It is left as a future work to examine post-Newtonian perturbations to triangular config- 
urations for general masses. The configuration may be non-equilateral or non-periodic. 
This work was supported in part (H.A.) by a Japanese Grant-in-Aid for Scientific Research 
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